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Abstract. 



> 

l/^ ' The distribution of capture zones formed during the nucleation and growth of point islands 

C3 ■ on a one-dimensional substrate during monomer deposition is considered for general critical 

island size i. A fragmentation theory approach yields the small and (for i = 0) large size 
asymptotics for the capture zone distribution (CZD) under the assumption of no neighbour- 
neighbour gap size correlation. These CZD asymptotic forms are different to those of the 
Generalised Wigner Surmise which has recently been proposed for island nucleation and 
growth models, and we discuss the reasons for the discrepancies. 



j^ . AMS classification scheme numbers: 45J05; 45M05; 78M35 



1. Introduction 



Island nucleation and growth during the early stages of thin film deposition has attracted 
much attention over the last few decades [28| HI ^7\. The islands that arise from the 
aggregation of the deposited material form the building blocks for nanostructure growth 
[211 [12] and subsequent film morphology [16], and so are of direct technological interest. The 
growth process is amenable to Monte Carlo (MC) simulation, which in turn have inspired 
many theoretical analyses which revisit meanfield rate equations (for example see references 



). The most striking aspect to arise from these models has been the identification of 
scale invariance in the growth process, exemplified in the island size distribution (ISD) [I]. 
It has been recognised that in order to understand the ISD, one cannot simply rely on mean 
field equations since they do not always yield the correct form for it [6] and in some cases do 
not even display the fundamental scaling property seen in simulation and experiment 



In order to better understand the ISD, we must consider the growth rate of an island that 
has nucleated on the substrate. It grows by capturing deposited monomers, which diffuse to 
it across the substrate. The growth rate of the island is therefore dependent on its capture 
zone, defined as the region of substrate closer to this island than to any other [TSl dSl [TT] : 
monomers deposited into the island's capture zone are more likely to diffuse to it than to any 
other competing island. This concept has been successfully used for the ISD found in growth 
on both two- [HI |6] and one-dimensional substrates [7J. The latter model is of particular 
interest, since it can be analysed in great detail, and it is the subject of the current paper. 

Recently there has been renewed interest in the capture zones and associated size distribution 
(CZD) where it has been proposed that the CZD follows the Generalised Wigner Surmise 
[2T1 122] . In particular, it has been suggested [21] that the CZD obeys 
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Here s = x/{x) is the capture zone size x scaled to its average (x) at any time t in the 
aggregation regime |2|. The aggregation regime is such that the island density is greater 
than the monomer density, so that further nucleation is a slow process compared to the 
island growth due to monomer capture. 

Despite the excellent visual comparisons between equation ([1]) and MC data taken from the 
literature [2T], the validity of the GWS has already been challenged. Shi et al. [25j studied 



i = 1 point island models in d = 1, 2, 3, 4 dimensions, finding that the CZD is more sharply 
peaked than the GWS suggests ior d = 1,4 and a better choice of /3 = 3 rather than (3 = 2 
for d = 2,3. Li et al. [H] also question the i = 1, d = 2 GWS form for the CZD, proposing 
their own form based on a sophisticated theory for capture zone evolution in two dimensions 

We also note that there is a confiict in the small and large size form of ([I]) and that of 
the theory presented in [7] for the one-dimensional (1-D) nucleation and growth of point 
islands. The work in [7] focused on the case with critical island size i = 1, where we impose 
the condition that z + 1 monomers must coincide in order to nucleate a stable, immobile 
island. The analysis in [7] utilises a fragmentation theory approach (as described below) 
which we consider to be physically reasonable. It is therefore of interest to ask whether a 
similar confiict exists for other critical island sizes. In this work we extend the 1-D point- 
island fragmentation approach to i = 0, 1, 2, 3, . . . and show that the small-size confiict with 
[2T] holds for all non-negative integers i. The case of i = 1 in this model has recently 
been considered in [13], where the form of equation ([1]) has been extended in light of MC 
simulation data. Here we retain a focus on the comparison between the original GWS of 
equation ([1]) and the fragmentation model for general i. 

In [7] a MC simulation was presented, where monomers are randomly deposited at rate F 
onto a line of points which represents the substrate. The monomers diffuse (with diffusion 
constant D) by nearest neighbour hops until they nucleate a new island, or are adsorbed by 
hopping onto an existing island. Since nucleation is a rare event in the aggregation regime, 
it is assumed [7] that the average monomer density ni{x) in the gap (length y) between two 
neighbouring islands obeys its saturated form 

ni(x) = Rx{y — x), < x < y, (3) 

F 
where we have set R = — . 

2D 

Furthermore, this density profile is then used to predict the probability of a nucleation 

occurring at position x in the gap, which we take to be proportional to n\^^(x). In adopting 

this form, it is assumed that the nucleation arises from the congregation of {i + 1) "mature" 

monomers that have explored their gap and lost memory of their deposition events. 

In this work we will adopt the same approach as in ^ and derive a gap-size evolution equation 
for i = 0, 1, 2, . . .. From this equation, we obtain information on similarity solutions which 
involve a reduced gap size distribution function 4>{y). We are then able to construct the 
CZD under the assumption used in [7] that due to the stochastic nature of the nucleation 
process there are no neighbour-neighbour gap size correlations in the system. Under this 



assumption, we can write 

r-2s 



P{s) = 2 r <i){y)<i){2s-y)dy. 
Jo 



(4) 



We shall then show that the small size form for P{s) given by equation (jl]) is different to 
that of equation ([1]). In particular, the small size power dependence from equation (jlj) will 
prove to be an odd power of s, whereas equation ([T]) always gives an even power for one- 
dimensional substrate. We also prove that the large- size forms in equations ([1]) and dl]) for 
the spontaneous nucleation model (i = 0) are in similar disagreement. 

2. The gap evolution equation 

We model the one-dimensional aggregation regime by regarding nucleated islands as points 
on a line. Our aim is to obtain an equation that describes the evolution of the gaps between 
adjacent islands. Let n{x,t) denote the number concentration of gaps of width x at time t. 
Since any new nucleation event that occurs in a 'parent' gap of width y will result in the 
creation of two 'daughter' gaps of widths x and y — x, it is clear that the evolution of the 
system of gap sizes can be interpreted as a fragmentation process. Processes of this type are 
usually modelled by the linear fragmentation equation 

d f°° 

— n{x,t) = -a{x)n{x,t)+ h{x\y)a{y)n{y,t)dy. (5) 

In the case of gap sizes, a{x) represents the rate at which gaps of width x fragment (due 
to nucleation) and h{x\y) describes the distribution of gaps of width x resulting from the 
fragmentation of a parent gap of width y > x. Since each new nucleation leads to the birth 
of only two daughter gaps, and since the total length of all the gaps is preserved, we require 
b to satisfy 

y rv 

b{x\y) dx = 2 and / xb{x\y)dx = y. (6) 

Jo 

To obtain appropriate functions a and b, as explained above, we follow the approach used 
by Blackman and Mulheran [7]. Therefore, we assume that, in the aggregation regime, the 
density of monomers in a gap of width y is given by ([3]) . 

By the reasoning given in the Introduction, the probability of a new nucleation occurring in 
a gap of width y may be taken as being proportional to 

a{y) = R+^ I x'+\y - xy+^ dx = R+^B{i + 2, i + 2) y^'+\ (7) 

^0 



where B{-,-) denotes the Beta function. Then, given that a nucleation event has caused a 
gap of width y to fragment, with fragmentation rate a{y), the probabihty that it will occur 
at a scaled position r = x/y in [ro, tq + dr] C [0, 1] is given by h{r)dr, where 

h{r) = R'+V+\l - ry+\ < r < 1. (8) 

Since x G [xq, Xq + dx] <^ x/y G [xo/y,XQ/y + dx/y] 4^ r E [ro,ro + dr], the probability that 
the nucleation occurs at x G [xq, xq + dx] is given by h{xo/y)dx/y. Therefore, we take 

b{x\y) = kh{x/y)/y, (9) 

where 

so that ([6]) is satisfied. This leads to 

b{x\y) = y-^'-^x'^\y - xy+^/B{i + 3, z + 2), 

and equation ([5]) therefore becomes 

d f 

— n{x, t)= R'+^ i -B{i + 2, i + 2) x^'+^n{x, t) 

+2^ x'+\y-xr+^n{y,t)dy\. (10) 

A simple re-scaling of the time variable then yields the gap-size evolution equation 



^ n{x, t) = -B{i + 2, i + 2) x^*+^n(x, t) + 2 / x'^'{y - xy+^n{y, t) dy. (11) 
ot Jx 

This equation is a particular case of the linear, homogeneous fragmentation equation 

d f°° 

—n{x, t) = -c^x'^nix, t) + c^ y''-^K{x/y)n{y, t) dy (12) 

for some cr > in which K{x/y), a homogeneous function of degree zero, determines 
the number of daughter "particles" of size x obtained when a parent "particle" of size 
y > X fragments. There is a considerable literature on equations of this type, and various 
mathematical techniques have been used in the analysis of f lT^ . For example, the theory of 
strongly continuous semigroups of operators can be applied to establish the existence and 
uniqueness of physically meaningful solutions; see, for example, [U Chapter 8]. Similarity 
solutions have also been investigated by a number of authors, including Ziff and McGrady 



[30] (for the case i = in equation ( ITTj) ). Cheng and Redner [9] and Treat [27]. As shown in 
[27] , similarity solutions can be written in the form 

,, , N*'^(t) , fN*(t)x\ ,^ , 

n*(x,t) = — Ai0[^_A_^j (13) 

where the reduced distribution is required to satisfy an integral equation and is normalised 
so that 

<P{y)dy= / y(t){y)dy= 1, 

^0 

and 

/•oo /'OO 

N*{t):= j n*{x,t)dx, V := / xn*{x,t)dx, 
Jo Jo 

are, respectively, the zeroth and first moments of n*. An explicit expression, involving 
the Meijer G-function, is derived in [271 Section 6] for the specific case when the daughter 
distribution function K in equation (fT2l) takes the form 



K{r)=r^{bo + b,r + --- + bprP), (14) 

p is a non- negative integer, bo, . . . ,bp eM., and < r < 1. We shall make use of the simple 
case p = I, discussed in [271 Section 7.2], in the next section. 

Asymptotic properties of the function have also been established for more general 
homogeneous functions K. In particular, in [9] and [271 Section 5], it is shown that 

0(y) = 0{y^) as y ^ 0, (15) 

provided that limr"'''^^ / sK(s)ds exists and is non- zero, and 
^^0 Jo 

(f)(y) = Oy^^^'>-^ expi-cy") as y ^ oo, (16) 

for some constant c > 0. 

The questions of existence and stability, in an appropriately defined sense, of similarity 
solutions to fragmentation equations of homogeneous type have also been addressed in [TO] . 

In the case of the gap fragmentation equation (fTTj) . (TT5l) and flT6|) lead directly to the following 
result. 

Theorem 1. The equation fll\) has a similarity solution of the form (Q^j where the gap size 
distribution satisfies 

6 



I. <P{y) = 0{f+^) asy^ 0; 

II. <p{y) = 0(?/~^exp(— ct/^*+^)) as y —7- oo for some c > 0. 

The task now is to understand, given the information we have for (f){y), the behaviour of the 
convolution (j4]). The situation for small s is straightforward for any i, but for large s we 
only have a result if z = 0. The derivation of the asymptotics is non-trivial and requires the 
use of a multi-dimensional version of Laplace's method. This is the content of the next two 
sections. 



3. CZD asymptotics 

The small s behaviour of P{s) for arbitrary i can be obtained immediately from part I of 
Theorem [1] and equation (jlj): 

Theorem 2. For the critical island size i > 0, we have 

P{s) = 0{s^'+^) ass^O, 
where P{s) is the 1-D CZD. 

Proof. We have, for small s and < y < 2s, 

0(y) = 0(y*+i), 0(2s -y) = 0((2s - yf^'), 
Hence (j){y)(j){2s — y) = Oi^s^'^^y^^^) and therefore 

P{s) = 2 / 0(y)0(2s -y)dy = 0{s'+^s'+^) = 0{s^'+^) 
Jo 

as s -)■ 0. D 



By Theorem [2l the exponent is always odd which differs from the GWS prediction that we 
should have P{s) = (9(s'^) where /3 = 2(i + 1) is always even when d = 1. Hence we are led 
to the conclusion that the GWS does not describe the behaviour of the scahng function P(s) 
as s — 7- for any z > if we accept that the gap evolution equation (ITTll . and the relation 
between (f){y) and the capture zone distribution P{s) given by (jl]), are correct. 

The next aim is to understand the asymptotic behaviour of P{s) in (j4]) for large s. For this, 
it would appear that an explicit and reasonably simple, formula for the reduced distribution 
(j){y) is required. In the case i = 0, the gap fragmentation equation (ITTj) becomes 

x{y-x)n{y,t)dy, 
7 



d , , 


x3 , r 


Qn{x^) = 


~n(.,t) + 2l 



which is the equation analysed by Ziff and McGrady in [30] . Note also, that in this case the 
homogeneous function K can be obtained from the general linear daughter distribution, 

K{r)=r^{bo + bir), (17) 

investigated by Treat in [271 Section 7.2], by setting 7 = l,6o = 12, 61 = —12. On applying 
[271 (7-6)], we deduce that (f){y) = 0(77), where rj = y^/fi^, 

_ r(2/3)r(7/3) _ 4 

^" r(4/3)r(2) -3^^^/^^ 



and 



Hence 



/^ r(i) Jo 

A f°° A /"°° 

= _7yV3 / z~^/'^ exp{-zr]) dz = —T]'^^^ / u~^^^ exp{-u) du. 

A /"°° 

Cf){y)= y^ U~^/^ exp{-u) du (18) 

^^ ^"t;-4/3e-^(^/'^)'cit;. (19) 



r(i)^2 A 

Note that Ziff and McGrady proposed a similar formula for but with the constant 
^^{k. 5.88) in the lower limit of integration in (I18p replaced by 6. That the correct choice is 
11^ follows from the fact that this leads to 

<P{y)dy= / y(t){y)dy= 1. 



This explicit representation of 0(|/) for the case i = enables us to obtain the following 
result. 

Theorem 3. If i = then 



P{s) = O fs-'/' exp (-2^)] 



as s —)■ oo. 



The next section is devoted to proving this theorem. 



4. Proof of Theorem [3] 

The arguments of Wong [29l Ch. VIII, Sections 8 and 11] allow us to obtain the following 
general theorem about two-dimensional Laplace integrals. 

Let 

J{\)= ff g{v,w)e-^f^''^'"^dvdw, (20) 



J J D 

where A is a large parameter, D C M^, and let f{v,w), g{v,w) be smooth functions on D. 
We have 

Theorem 4. // (t>o, wq) is the global minimum of f{v, w) on D and is a critical point of the 
third kind, then as X^oo, J(A) = C(A-2g-A/(t,o,«^o)), 

We remind the reader that {vq^Wq) is a critical point of the third kind of f{v,w) if it 
is an extremum point of / on D belonging to the boundary of D through which pass two 
intersecting tangent lines, neither of which coincides with the tangent to the level curve 

ip{v,w) = ip{vo,Wo) at {vo,Wo). 

Our aim is to understand the asymptotics of (jl]) as s — )■ cxd in the case of z = using the 
explicit form f lT9|) . 

If we let u = y^v/fi^ in ( ITSl) . substitute the resulting expression for (j){y) into (jl]) and then 
put y = 2sz, we have 

/i4r(2/3)^ Jl Jl Jo 

Our strategy is to use Laplace's method to evaluate the inner integral, which will give us 
a two-dimensional Laplace integral on D = (1,cxd) x (1,oo), which will be attacked using 
Theorem HI 

For any {v,w) G D, set S{z) = z^v + (1 — zYw. li v = w, S{z) has a unique minimum at 
z+ = 1/2. li V j^ w, critical points of S{z) satisfy {v — w)z^ + 2wz — tf = 0, so that 

—W ± y/vW 
Z± = . 

V — w 
It is easy to check that 5"'(z+) = 6y/vw > 0, so that the minimum of S{z) is obtained at 



z+ = ^L_ = y e 0, 1 . 22 

V — w Jw + Jv 



8s^ 
Hence by Laplace's formula, setting A = — -, we have, as A — )■ cxd, 

/i3 



Ji(A) := ^' zil - z) expi-XSiz)) dz ~ z^l - ^+)e-^^(^+) J ^-^^. 



Note that 

wv 



which reaches its niininium value 1/4 on D at the corner point (1, 1). 
Thus, the inner integral in (1211) satisfies, as s — )■ oo. 



^^{—r]^{^r] P(^' ^) exp 



where we have put 



Piv^^) = \l-gJi^^+i^-^+)- 



This means by (1211) that 



^<'> ^ ° (^"^ r /°° "<■''• '"><■'"">"" ^^p {-'iu^rvw) *""") ■ •''' 



and now invoking Theorem H] we have, as s — )■ cxd, 

o 3\ -2^ 



i.e. P{s) = O i s~^/^e~^'* ^'^ ) as required. 



5. Conclusions 

In this paper we have shown that the Blackman and Mulheran model |7j for the nucleation 
and growth of point islands on a one-dimensional (1-D) substrate, based on a fragmentation 
equation approach, yields a capture zone distribution (CZD) that is different to the 
Generalised Wigner Surmise (GWS) in equation ([1]) [21]. The asymptotics of the CZD 
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are measurable in MC simulation, allowing these theories to be tested, and we will present 
our analyses of simulations elsewhere [20] • Here we will conclude this current work with a 
discussion of how the two theoretical approaches differ and how they might be reconciled. 

It is interesting to note that the justification for the relationship between the parameter /3 
in equation ([1]) and the critical island size, as given in equation ([2]), is based on the same 
physical model analysed in this paper. In [21], the island nucleation rate is discussed in 
terms of the monomer density n, and the probability of {i + 1) monomers coinciding is used 
to give the nucleation rate as ~ ra*"*"^. This of course is exactly the physical basis we have 
used, so it is worthwhile considering why [2T] end up with different power-law behaviour for 
small capture zone sizes. 



Let us summarise the phenomenological arguments used in [2T]. Firstly the authors consider 
the nucleation rate within capture zones with size s, which have a relative density P{s), 
using the mean field value of the monomer density n. They then reconsider the nucleation 
rate within a small zone of size s using the locally averaged monomer density. In one space 
dimension, they argue that this density ~ ns"^. Equating the two rates, P{s)n^~^^ ~ (ns^)*+^ 
yields their law for /3 in 1-D nucleation so that P{s) ~ s^*+^ for small s. We note that 
recently the originators of equation ([1]) have considered alternative forms to the Gaussian 
tail following evidence from their own MC simulations with i = 1. 

From our perspective, we see a number of problems with this argument. Firstly, the method 
of constructing the two alternative nucleation rates appears internally inconsistent, using 
both the mean field monomer density n and a local approximation. Usually this approach 
would involve taking averages of the latter to reach consistency with the former [3l [7], 
however this was not done in pi]. The second criticism, more pertinent to this paper, is 
that the nucleation rate using the locally averaged monomer density within a capture zone 
has not been justified. Indeed, in [22], the authors revise their argument for two dimensional 
substrates and use a spatially dependent monomer density n{r) within a capture zone to 
derive the nucleation rate for the zone. This involves integrating n*"'"^(r) over the zone, 
reminiscent of the approach used elsewhere for the two dimensional substrate [15] , resulting 
in a different relationship (3{i) which better fits Monte Carlo datajj 

Presumably, this type of modification can be taken forward to the 1-D case discussed here. 
Following the same line of argument, we might take the local monomer density in a small 
capture zone size s as that from two gaps of size s/2, using the form of equation ([3]) above. If 
we do this, we will end up with the formula I3{i) = 22 + 3, thereby agreeing with the small-size 

I In ;22], the authors use n{r) ^ R^ — r^ , Ri < r < R, with Ri and R being the island and circular capture 
zone radii. It is not clear what motivated the choice of boundary conditions for this form; we would have 
expected the monomer density to be zero at the island edge, rising to a maximum at the capture zone radius 
so that its gradient is zero at this boundary [15]. 
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CZD behaviour we have derived in this work. Although this does go some way to reconcihng 
the two approaches, we consider that this hne of reasoning is heuristic at best, and much 
prefer the more rigorous methods adopted in our work presented here in this paper. 

The second component of the GWS in equation ([1]) is the large-size behaviour of the CZD. 
Here too we see a conflict with the fragmentation theory approach we adopt, again despite 
the same physical basis for the models. In [21] , the authors are motivated by the idea of a 
large fluctuating capture zone being constrained by its neighbours, rather than any attempt 
to consider in detail how the nucleation process impacts on the large size CZD behaviour. 
In this work we have been able to prove that for z = 0, the large size dependence of the 
CZD following equation (jl]) mirrors that of the corresponding gap distribution 0(|/). In 
particular, we have P{s) ~ exp{—2s^/^^) rather than the Gaussian tail of equation ([T]). 
Furthermore, whilst it remains unproved, we can conjecture that the CZD for larger values 
of i will similarly follow the asymptotics of the corresponding GSD and follow exp(— qs^*"^^) 
for some constants q in stark contrast to the Gaussian tail. 

In conclusion, we have shown that the GWS of |21] is in conflict with the asymptotic solutions 
to the fragmentation theory analysis of point island nucleation and growth in one dimension. 
Given that the physical basis of the two approaches is the same, we believe that these 
differences show important failings of the GWS. Of course, this is not to claim that the 
Blackman and Mulheran model, embodied in equations ([3]) and @ above, has been proven 
to be correct; confrontation with simulation and experiment will ultimately arbitrate between 
these theories. Along these lines, we note the recent analysis of the case of i = 1 in [13] , and 
our own detailed comparisons with extensive simulation data for i = 0, 1,2,3 presented in 
101. 
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